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We show that superconductors with broken time-reversal symmetry have very specific magnetic 
and electric responses to inhomogeneous heating. A local heating of such superconductors induces 
a magnetic field with a profile that is sensitive to the presence of domain walls and crystalline 
anisotropy of superconducting states. A nonstationary heating process produces an electric field 
and a charge imbalance in different bands. These effects can be measured and used to distinguish 
s + is and s + id superconducting states in the candidate materials such as Bai_ a; K a; Fe 2 As 2 . 
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In many recently discovered superconducting materi¬ 
als, the pairing of electrons is supposed to take place in 
several sheets of a Fermi surface formed by overlapping 
electronic bands [1-6]. Of special interest are the states 
where the difference of gap’s phases in the bands is nei¬ 
ther 0 or 7T [7-18]. Indeed, in addition to the breakdown 
of usual U(l) gauge symmetry, such superconducting 
states are characterized by an extra broken time-reversal 
symmetry (BTRS) that has numerous interesting phys¬ 
ical consequences, many of which are not yet explored. 
Iron-based superconductors [3] are among the most com¬ 
monly accepted candidates for the observation of a BTRS 
state originating from the multiband character of super¬ 
conductivity and several competing pairing channels. 

Experimental data suggest that in the hole-doped 122 
compounds Bai_ :c K :c Fe2As2 the symmetry of the super¬ 
conducting state can change depending on the doping 
level x. At moderate doping x ~ 0.4 various measure¬ 
ments including neutron scattering [19], thermal con¬ 
ductivity [20], and angle-resolved photoemission spec¬ 
troscopy (ARPES) [21-23] are consistent with the hy¬ 
pothesis of the s± state where the superconducting gap 
changes sign between electron and hole pockets. On the 
other hand, the symmetry of the superconducting state 
at heavy doping x 1 is not so clear regarding the ques¬ 
tion of whether the d channel dominates or if the gap 
retains s± symmetry changing sign between the inner 
hole bands at the T point [24, 25]. Indeed, there are ev¬ 
idences that d-wave pairing channel dominates [26-29], 
while other ARPES data were interpreted in favor of an 
s-wave symmetry [30, 31]. 

Which of these two possibilities is realized at heavy 
doping depends on the fine balance of the pairing interac¬ 
tions in different channels. However, both cases strongly 
suggest the existence of an intermediate superconducting 
state that breaks time-reversal symmetry at a certain 
range of the doping level x. Two alternative scenarios 


* Electronic address: garaud.phys@gmail.com 



Figure 1: Variation of the interband phase difference 0 12 
in BTRS three-band superconductors induced by a hot spot 
created, e.g., by a laser pulse. The phase difference variation 
induced by temperature gradients around hot spots in the 
case of an s + is state (left) preserves Ca symmetry, while it 
has fourfold structure for an s + id state. The value of Oo in 
the s + is case is 0.18 while for s + id it is smaller: 60 — 0.05. 


have been considered: namely, s + id and s + is symme¬ 
tries [8, 9, 17, 30, 31]. The s + id state is anisotropic, 
as it breaks C4 crystalline symmetry, while the s + is 
state is qualitatively different as C4 symmetry is pre¬ 
served [17]. Note that the s + id state is qualitatively 
different from the (time-reversal-preserving) s + d state 
which earlier attracted interest in the context of high- 
temperature cuprate superconductors (see, e.g., [32-35]). 
It also contrasts with the d + id state that violates both 
parity and time-reversal symmetries [7, 36] . 

To this day, no experimental proof of s + is or s + id 
BTRS states has been reported. Indeed, probing the rel¬ 
ative phases between components of the order parameter 
in different bands is a challenging task. For example, 
the s + is state does not break point group symmetries 
and therefore it is not associated with intrinsic Cooper 
pair angular momentum. Hence, it cannot produce a lo¬ 
cal magnetic field and is invisible for conventional meth¬ 
ods like muon spin relaxation and polar Kerr effect mea¬ 
surements that were used to search for BTRS p + ip 
superconducting state in, e.g., the Sr2Ru04 compound 
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[37]. Proposals for indirect observation of a BTRS sig¬ 
nature in pnictides, with various limitations, have been 
recently voiced. These include, for example, investiga¬ 
tion of the spectrum of collective modes which includes 
massless [14] and mixed phase-density [15, 17, 38, 39] ex¬ 
citations. It was also proposed to consider exotic topo¬ 
logical excitations in the form of skyrmions and domain 
walls [40-42] , an unconventional vortex viscosity mecha¬ 
nism [43], vortex clustering [15], and exotic reentrant and 
precursor phases induced by fluctuations [44-47]. Spon¬ 
taneous currents were predicted to exist near impurities 
in anisotropic superconducting states [ 8 , 18] or in sam¬ 
ples subjected to strain [18]. The latter proposal actually 
involves symmetry change of s+is states and relies on the 
presence of disorder, which usually has an uncontrollable 
distribution. 

In this Letter, we discuss an experimental set-up based 
on a local heating, that allows the direct observation of 
BTRS states in a controllable way. This is illustrated 
in Fig. 1, where local heating induces a local variation 
of relative phases that are further shown to yield elec¬ 
tromagnetic excitations. The key idea is based on the 
recent proposal of an unconventional thermoelectric ef¬ 
fect in BTRS multiband superconductors [48]. There, a 
temperature gradient generates phase gradients of con¬ 
densate components, due to the generically temperature- 
dependent interband phase differences 0ki(T) := Ok — 0{ 
(where k,i are band indices). It results that the local 
heating generates spontaneous magnetic fields and charge 
imbalance distributions. These thermoelectric responses 
are drastically different from their counterparts in con¬ 
ventional superconductors [49]. As discussed below, the 
fields created by local heating have opposite directions 
in two degenerate superconducting states [i.e. s + is(d) 
and s — is(d ) ones]. They are measurable by conven¬ 
tional techniques (e.g. by SQUID) and therefore scans 
of the surface can be used to diagnose the structure of 
order parameter and interband phase differences, to de¬ 
tect pinned domain walls or broken crystalline symmetry 
states in either s + is or s + id superconductors. 

We consider a minimal three-band microscopic model 
which has been suggested to describe the BTRS super¬ 
conducting state in the hole-doped 122 iron-pnictide com¬ 
pounds [10, 17, 39] with three distinct superconducting 
gaps Ai 5 2,3 in different bands. The pairing which leads 
to the BTRS state is dominated by the competition of 
different interband repulsion channels described by the 
following coupling matrix: 

/O V x\ 

9 = \ V 0 A . ( 1 ) 

\ A A 0 / 

Here, we assume for simplicity that the density of states 
vo is the same in all superconducting bands. This model 
has been suggested [10, 17, 39] in order to describe transi¬ 
tions between s/s± and s+is states when tuning parame¬ 
ters 77 ,A and temperature. The dimensionless coefficients 
77 and A describe different pairing channels, whether it 


is an s + is or an s + id state. In the former case Ai 5 2 
correspond to the gaps at hole Fermi surfaces and A 3 is 
the gap at the electron pockets, so that u^h = ^oV and 
u e h = vo\ are, respectively, the hole-hole and electron- 
hole interactions [17, 39]. The same model (1) can be 
used to describe the s + id states but there, Aq 2 describe 
gaps in electron pockets and A 3 is the gap at the hole 
Fermi surface, so that u e h = voX and u ee = vof] are 
electron-hole and electron-electron interactions respec¬ 
tively. 

To study magnetic and electric responses of both s + is 
and s + id states, we use a Time-Dependent Ginzburg- 
Landau (TDGL) approach [50, 51] generalized to a multi¬ 
band case (see derivation in Appendix A and B). The 
dimensionless TDGL equations read (see details in [67]): 

dXF 

(d t + 2ietp)ip k = -—, V x B - a n E = j s , ( 2 ) 

°Vk 

where ip is the electrostatic potential, a n is the normal 
state conductivity, and j s = —dTjdA is the supercon¬ 
ducting current. Near the critical temperature the energy 
relaxation is determined by the phonon scattering which 
yields the relaxation time scale to = nh/(8T c ) ^ 1 ps 
provided T c ~ 1 meV, which is about 10 K. 

Note that multiband superconductors are described by 
several components ipk which do not necessarily coincide 
with the gap functions A^ in different bands (see, e.g., [17, 
39]). For example, since the coupling matrix (1) has only 
two positive eigenvalues, the relevant Ginzburg-Landau 
(GL) theory reduces to a two-component one (see details 
in Appendix A), which in dimensionless units reads as 

T = + 5Z ( k 3 l n ^I 2 + a 3 l^iI 2 + y l^jI 4 ) ( 3a ) 

3 = 1 ' 

+ fci 2 , a (n;^n a ^2 + c. c .) (3b) 

+ 7lV’i| 2 |V’2 | 2 + ^(V’i 2 V>2 +C.C.) , (3c) 

with II = V — 2 ieA. The components Vh, 2 are deter¬ 
mined by a superposition of the different gap functions 
A i. All coefficients of the model (3) are consistently de¬ 
termined from the microscopic coupling matrix (see de¬ 
tails in Appendix A), and the temperature dependence 
is given by the coefficients: 

aq = —2(Go — G\ + t) (4a) 

c ^2 = —(2x 2 + l)(Go — G 2 T t) (4b) 

where r = (1 - T/T c ), x = (rj - +8A 2 )/(4A), Gi )2 

are the positive eigenvalues of the inverse coupling matrix 
vq g~ x and Go = min(Gi, G 2 ). The general GL functional 
(3) derived from the three-band microscopic model has 
S > 0. Hence it favors BTRS with +tt/2 phase differences 
between 'ipi and ip 2 order parameter components, which 
describes both the s + is and s + id states depending on 
the structure of mixed gradient terms (3b). They are 



ki2,x = ki 2 , y fc> r the s + is state and ki2,x = —&i 2 ,y for 
the crystalline (^-symmetry breaking s + id state. 

As a consequence of the discrete degeneracy due to 
BTRS, the model (3c) allows domain walls (DW) inter¬ 
polating between regions with different relative phases. 
The direct observation of DWs in s + is(d ) states is chal¬ 
lenging. Unlike DWs in p + ip superconductors, they 
do not generate a spontaneous magnetic field. However, 
by our general argument below, DWs should provide a 
controlled magnetic response in the presence of relative- 
density perturbations that can be induced by a local 
heating. 

To investigate the response to spatial modulations of 
the components of the order parameter induced by a local 
heat source, the fields ^ 1,2 and A are discretized using a 
finite-element framework [54] (see discussion of numerical 
methods in Appendix C). To model the local heating, the 
temperature profile is found by solving the (stationary) 
heat equation for a heat source at temperature T s , while 
boundaries are kept at Tq = 0.7T C . Once the temperature 
profile is found, the coefficient a*. in (3) varies in space 
and the TDGL equations ( 2 ) are evolved for At = 80 
(in units to defined above). The temperature of the heat 
source is then modified to T', and the TDGL equations 
are further evolved for the new temperature profile for 
a period At. The temperature of the source is initially 
set to To, sequentially ramped up to 0.95T C , and then 
ramped down back to Tq. In our simulations, we chose 
the dimensionless conductivity cr n = 0.1 and the coupling 
constant e = 0.113. The coefficients in GL functional (3) 
are determined using the microscopic coupling matrix ( 1 ) 
with coefficients 77 = 5 and A = 4.5 [ 68 ]. 

As shown in Fig. 2 , according to the simulations, when 
the source heats up a domain wall, it induces a multipo¬ 
lar magnetic field with zero net flux. In the case of a 
superconductor with s + is symmetry, it shows a dipo¬ 
lar structure, while it is differently distributed for s + id. 
On the other hand, when the heat source is focused on 
the uniform s + is state it shows no magnetic response, 
while a fourfold magnetic field is induced in the s + id 
case as a result of the explicit breakdown of the C 4 sym¬ 
metry. Here, spatial variations are normalized to the 
penetration depth A l and the amplitudes of the induced 
magnetic field to the second critical field H C 2 defined by 
the GL functional (3). Provided the typical values of 
H c 2 ^ 10T in 122 pnictides [56], one can see that the 
magnetic response can be detected with high accuracy 
by conventional local probes of static magnetic field such 
as scanning SQUID or Hall probe microscopy. 

The physical origin of the spontaneous magnetic re¬ 
sponse follows from that the total current is the sum of 
partial currents in each of the N bands j = J2^ =1 jk 
and therefore can generated by the gradients of relative 
phases [48]. Since jk = (V0& — 27 rA/ 4 > 0 )c 4 > 0 /( 87 r 2 A^) 
the London expression for the magnetic field in multi- 
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Figure 2: The magnetic response that originates in local 
heating of the sample. Panels (a) and (c), respectively, show 
the response when a local hot spot heats an area of a sample 
which contains a domain wall. Here, we display two cases of 
the domain walls for two BTRS states: s+is and s+id. Panels 
(b) and (d) show the response to the local hot spot in the 
case of the homogeneous BTRS state for the s + is and s + id 
superconducting state. The color plot shows the magnitude of 
the out-of-plane induced magnetic field B z (magnitudes differ 
in different panels), while arrows indicate the orientation of 
supercurrents. The dotted line indicates the presence of the 
domain wall and the inhomogeneous temperature profile is 
induced by the ringlike heat source shown by the small circles. 
Length scales are given in terms of London penetration length 
and calculated values of the GL parameters are given in [52, 
55] and Appendix A-C. 

band superconductors is modified as follows 

B = —■ v x (A \j) + JsL £ V x { lki V0 ki ) (5) 

k>i 

where A k are coefficients characterizing contribution 
of each band to the Meissner screening, A l = 

1 / \j^k A& 2 the London penetration depth and 

+ki( r ) = ^i(^fc 2 ( r ) — Ar 2 (r)). In contrast to London’s 
magnetostatics, Eq. (5) shows that the magnetic field 
features an additional contribution when relative density 
gradients V 7 ^(r) are noncollinear with that of relative- 
phase gradients . Such gradients generically appear 
in BTRS states if a domain wall-containing superconduc¬ 
tor is exposed to a local heat source. The second term in 
(5) can be nonzero even in the absence of domain walls 
due to direction dependent tensor coefficients 7 /^( 1 *) in 
anisotropic s + id states. 

Domain walls can be created by quenching the sample 
and stabilized by pinning or artificial geometric barriers 
[42, 57]. Yet it is also important to obtain the evidence 
of isotropic s + is states for homogeneous superconduct¬ 
ing states. Below, we show that this can be done by 
considering the nonequilibrium electric responses gener¬ 
ated by nonstationary heating when the local tempera¬ 
ture evolves recovering from the initial hot spot created, 
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Figure 3: Electric response due to the nonstationary heating 
of the sample. Panel (e) shows the time evolution of the 
source’s temperature. Panels (a)-(d) correspond to the cases 
in Fig. 2 and the color plot shows the voltage Vn generated by 
the charge imbalance that can be picked by a normal detector, 
while the arrows correspond to the electric field. Panel (e) 
also displays <f> integrated over the whole sample, of an s + is 
superconductor without DW (b). The red dot on panel (e) 
denotes the “position” of panels (a)-(d) in the time series. 


e.g., by a laser pulse [49]. An unusual electric response 
can be seen when combining Eq.(5) to Faraday’s law. 
The electric field E = —c~ 1 d t A — Vy? can be rewritten 
as 

E =! 

k>i 

Here, <f> = ^2 k {(p + hO k /2e)/N is a gauge invariant po¬ 
tential field, determined by the sum of chemical potential 
differences between quasiparticles fi q = etp and conden- 

sates in each band /q> = —hOk/ 2. Each of the partial 

potential differences = [fi q — is proportional 

to charge imbalance in the k -th band Q k = 2e 2 uo^ k ^ 
[58-60]. In multicomponent systems the charge imbal¬ 
ance can generated by variations of interband phase dif¬ 
ferences in space and time. The physics behind this pro¬ 
cess is a nonequilibrium redistribution of Cooper pairs 
between different bands which initially creates partial 
charge imbalances Q* k . This mechanism leads to the un¬ 
conventional electric response of BTRS superconductors 
to a nonstationary local heating. It can be measured 
with potential probe techniques that were employed to 
study the imbalance between quasiparticles and conden¬ 
sate subsystems in conventional superconductors [61, 62]. 

Fig. 3 shows such an electric response to a nonstation¬ 
ary heating of the superconducting sample. This multi¬ 


component electrodynamic phenomenon can be used to 
detect BTRS states through charge imbalance generation 
in response to nonstationary heating. As shown in Fig. 3, 
the charge imbalance shows nontrivial pattern that is dif¬ 
ferent for s + is and s + id states. The total charge im¬ 
balance (<F) in uniform s + id is zero as a result of the 
fourfold symmetric structure due to broken C 4 symme¬ 
try. On the other hand, as shown in Fig. 3(e), (<f>) 7 ^ 0 in 
uniform s-j-is. Note that the total imbalance in the case 
of domain walls Fig. 3(a) and Fig. 3(c) is zero because 
the heat source is centered at the domain wall. A shifted 
source from the DW center will not be symmetric and 
thus its average will not vanish. 

The generated charge imbalance can be measured with 
the standard technique using the normal metal and su¬ 
perconducting potential probes [61, 62]. The magni¬ 
tude of voltage Vjv induced in the normal detector is 
related to the dimensionless signal shown in Fig. 3 as 
Vjv = h<&/(toe), where h/(t oe) « 4 mV. The overall 
magnitude of the voltage signal is thus expected at the 
order of the fiV. The electric field and charge imbal¬ 
ance depend on the dynamics of the temperature profile 
variation, while B depends on the temperature profile it¬ 
self. As a result, the sign of the induced electric field and 
charge imbalance changes when ramping down the heat- 
source temperature, while it does not for the magnetic 
field (see animations [63] and their description in Ap¬ 
pendix D). Note that similarly to the spontaneous mag¬ 
netic field, the induced charge imbalances are sensitive to 
BTRS: degenerate s + is and s — is states produce oppo¬ 
site electric fields and charge imbalances in response to 
the same heating protocol. It allows to discriminate be¬ 
tween the usual thermoelectric occurring in conventional 
superconductors and the unconventional one being a spe¬ 
cific signature of BTRS states. 

To conclude, we demonstrated possible direct mani¬ 
festations of BTRS states in experimentally observable 
electric and magnetic responses to nonuniform and non¬ 
stationary heating. The signs of the generically induced 
magnetic field and charge imbalance distributions are op¬ 
posite in degenerate BTRS states [i.e. in s + is(d) and 
s — is(d)]. These specific thermoelectric behaviors were 
also shown to reveal the presence of domain walls between 
s + is{d)/s — is{d) states. Moreover, the demonstrated 
crucial dependence of thermomagnetic and charge im¬ 
balance responses on crystalline anisotropy provides an 
experimental tool to distinguish between isotropic s + is 
and C4 symmetry-breaking s + id states that are partic¬ 
ularly interesting for pnictides. 
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Appendix A: Ginzburg-Landau expansion for the 
s + is and s + id superconductors 

We consider here two alternative patterns of super¬ 
conducting coupling which both result on BTRS state 
but will be shown to yield qualitatively different phys¬ 
ical properties. In the first scheme shown in Fig. 4(a) 
the dominating pairing channels are the interband repul¬ 
sion between electron and hole bands, as well as between 
two hole pockets at T. There, the order parameter is the 
same in both electron pockets so that the crystalline C 4 
symmetry is not broken and thus corresponds to an s 
state. Instead the second alternative shown in Fig. 4(b) 
is that when the strongest interactions are the repulsions 
between hole and electron bands and between two elec¬ 
tron pockets. Such interaction favours order parameter 
sign change between electron pockets resulting in a C 4 
symmetry breaking d-wave state. 

To derive the Ginzburg-Landau expansion that is used 
in our simulations, we consider the microscopic model 
of clean superconductor with three overlapping bands at 
the Fermi level. Within quasiclassical approximation the 
band parameters characterizing the two different cylin¬ 
drical sheets of the Fermi surface are the Fermi velocities 
vP and the partial densities of states (DOS) 1 z,-, labelled 
by the band indices j = 1,2,3. To describe the two 
possible alternatives of BTRS states, namely s + is and 
s + id symmetries, we consider two three-band models 
schematically shown in Fig. 4. The Eilenberger equa¬ 
tions for quasiclassical propagators take the form 

hv^Ilfj+ 2u n fj-2A j g j =0, (A.l) 

hv^n*f+ - 2 L 0 n f+ + 2A*g 3 = 0, 

where II = V — 27riA/4>o, A is the vector potential, 
vP is the Fermi velocity. The quasiclassical Green’s 
functions in each band obey normalization condition 
g] + fjfj~ = 1 - The self-consistency equations for the 
gaps and electric current are 

Ai(p,r) = 2nT ^2 ^ij(p,P')fj(p,r,u n ) (A.2) 

n,p',j 

j(r) = 2 tt eTu Y v^lm gj(p, r,u n ) (A.3) 

n,P,j 




Figure 4: (Color online) - Schematic view of the band struc¬ 
ture in hole-doped iron pnictide compound Bai_ cc K cc Fe 2 As 2 . 
It consists of two hole pockets at T shown by open circles 
and two electron pockets at (0; n) and (77; 0) shown by filled 
ellipses. In panel (a), the s + is state is favored by supercon¬ 
ducting coupling dominated by the interband repulsion be¬ 
tween electron and hole Fermi surfaces Uhe , as well as between 
the two hole pockets Uhh • In panel (b), there is a possibility 
of s + id state due to the repulsion between electron and hole 
Fermi surfaces Uhe as well as between electron pockets u ee . 


the expansion by powers of the gap functions amplitudes 
A j and their gradients: 

fj(p,r,u n )= (A.4) 

UJ n 2c 2 c 0^ ' 4u 

and fj- (p, r, uo n ) = /*(—p, r, u n ). Then, for the summa¬ 
tion over Matsubara frequencies, we get 

N d 

27 tT ^ ^ uj n 1 = Go + t (A.5) 

n =0 

where r = (1 — T/T c ). We normalize gaps by T c /y/p, 
where 




7rT c 3 _ 7C(3) 


n=0 




87r 2 


0 . 1 , 


(A.6) 


substitute (A.4) into the self-consistency Eqs.(A.2) and 
get the system of GL equations 

[(Go + r — A" 1 ) A] . = -A'^ ) n a n 6 A i + |A i | 2 A i , (A.7) 


where gj = sign(cj n )y 1 — fjfj~ and v is the density of 

states, p parameters run over the corresponding Fermi 
surfaces and \j is the coupling potential matrix. For 
simplicity we will consider further isotropic pairing states 
so that A ij(p,p') = const on each of the Fermi surfaces. 
However in electron pockets we keep the anisotropy of 
Fermi velocities in Eq.(A.l). We neglect the anisotropy 
of hole bands which is a well-justified assumption [64]. 

The derivation of the GL functional from the micro¬ 
scopic equations formally follows the standard scheme. 
First we find the solutions of Eqs.(A.l) in the form of 


Here A = (Ai,A 2 ,As) t , and the anisotropy tensor is 

= h2 P( v Pa v Fb)/^Tc where the average is taken over 
the j- th Fermi surface and a, b stand for the x, y coordi¬ 
nates. The current is given by 

j(r) = Y , Im V/GIIA,;. (A.8) 

^ i 

In the following, we consider coupling matrix g = no A 
describing the case of an interband dominated pairing 
with repulsion. We assume that the density of states is 
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the same in all bands ^1,2,3 = v 0 and parametrize dif¬ 
ferent pairing interactions with two dimensionless coeffi¬ 
cients 77 and v as follows 


/ 0 T] X\ 

A = - 7? 0 A . (A. 9 ) 

\A A 0 j 


tions 


aiipi +b lj \ip j \ 2 ipi +bjiplipl = (A. 13 a) 

(KW + k( 2 J )+ x(Ki 2 J - K$) n ^ 2 
a 2 il > 2 + b 2 j\ipj 1 2 V>2 + bjtp^i = (A. 13 b) 



Here the coefficients describe different pairing channels 
in s± + is and s + id states. In the former case Ai^ 
correspond to the gaps at hole Fermi surfaces and A3 
is the gap at electron pockets so that u k h = v orj and 
u e h = vo \ are the hole-hole and electron-hole interac¬ 
tions correspondingly. In contrast to describe the s + id 
state we use the same model (A. 9 ) but assume that Ay2 
describe gaps in electron pockets and A 3 is the gap at the 
hole Fermi surface so that u ek = vod and u ee = z^A to 
be electron-hole and electron-electron interactions corre¬ 
spondingly. 

Neglecting the r.h.s. in (A. 7 ) we get the lin¬ 
ear equation which determines the critical tempera¬ 
ture Go = min(Gi,G2), where G 1 = 1 /77 and G2 = 

(v + vV + 8A 2 ) / 4 A 2 are the positive eigenvalues of the 
matrix 


A- x = 


1 

2 X 2 r] 


A 2 -A 2 —A77 \ 
—A 2 A 2 —A rj J 
—A r] —Ary t ] 2 J 


(A.10) 


The coupling matrix A 1 has only two positive eigenval¬ 
ues Gi ? 2 whose eigenvectors are Ai = (—1,1,0) T and 
A2 = (x,x, 1 ) T with x = (rj — yjif + 8 A 2 )/ 4 A. Since 
only the fields corresponding to the positive eigenvalues 
can nucleate, the GL theory (A. 7 ) has to be reduced to 
a two-component one. To implement this reduction we 
represent the general order parameter in terms of the 
superposition 


The parameters of the left hand side of the Ginzburg- 
Landau equations (A. 13) are expressed, in terms of the 
coefficients of the coupling matrix (A.9) as 

aj = -\Aj\ 2 (Go — Gj+r), (A.14a) 

with | Ai | 2 = 2 and |A 2 | 2 = 2x 2 + 1 (A.14b) 

611 = 2 , 622 = (2x 4 + 1) and b k := b kk (A.14c) 

612 = 4x 2 , bj = 2x 2 . (A.14d) 

The system (A. 13) is quite general and describes both 
s-\-is and s + id states. The difference between these two 
cases is determined by the symmetry of gradient terms. 
In s + is state the G 4 symmetry requires that K X x = 
Kyy = K^\ Then the general Eqs. (A.13) are simplified 
as follows 

aiVh + 6ij|Vh| 2 Vh + bj^l = k^U 2 ^ 

CL2^2 + b 2 j\^j\ 2 ^2 + = k 2 jli 2 'ip j (A.15) 

where 

k lt =KW +KW 
k 22 = [(K (1) +K^)x 2 + K^} 

k 12 = x[KW -K(% (A.16) 

In s + id state the G 4 symmetry is broken so that 
K^ x 2 ^ = Kyy^ = K^’ 2 \ where the bands 1,2 corre¬ 
spond to the electron pockets. The hole band 3 is con¬ 
sidered to be G 4 symmetric so that K^ S J = Kyy = 

Then we obtain from Eqs. (A. 13) 


A = Vh Ai + ^2 A 2 , (A.11) 

so that 

(Ai, A 2 , A 3 ) = (xip 2 - i)i,x^ 2 + Vh, ^ 2 ) • (A.12) 

where ijji ^ 2 are the complex order parameter fields that 
have different interpretation depending on the system 
considered. In the s + is case 7/q and 7/^2 are the order 
parameter of s± pairing channels between two concen¬ 
tric hole surfaces and between hole and electron surfaces 
correspondingly. For the system with s + id symmetry, 
7/q is the order parameter of the d wave channel in elec¬ 
tron pockets and ^ 2 is the order parameter of s± pairing 
between electron and hole surfaces. 

Now, substituting the ansatz (A.ll) into the system of 
Ginzburg-Landau equations (A. 7 ) we obtain, after pro¬ 
jection on the vectors Ai^, the system of two GL equa¬ 


ai'ipi + hj | ipj 1 2 V>i + bjiplip% = (A.17a) 

knlf'ip! + k 12 H xy ip 2 

a2il>2 + b 2j | ipj\ 2 i> 2 + bjip^l = (A.17b) 

&22IIV2 + h 2 H-xylpl 

where the coefficients k ak are given by the same Eq. 
(A. 16) as in s + is state and U xy = U 2 — U 2 . 

The general free energy functional which gives both 
the GL equations (A. 15,A. 17) is thus given by F = ^ + 
BqF s where 

Fs = E { k n inV’il 2 + <*■1 2 + f IV’il 4 } 

3 = 1 

+ k 12 , a ((n a ^)*n Q ^ 2 + c.c.) + tIVTIVT 

+ ^ (^i 2 '*/ , 2 + c.c.) . (A.18) 
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We denote Bq = T c ^/z / 0 /p which is of the order the ther¬ 
modynamical critical field at low temperatures [53]. An 
extra factor of p here comes from our normalization of 
gaps. 

In the case of s + is symmetry, ki2,x = & 12 , y = &12 , 
while for s + id symmetry the mixed gradients satisfy 
k\2,x = ~k\ 2 , y = kw The other coefficients are given by 
= a k , Pk = b k k, 7 = h2 and S = bj. 


Appendix B: Time-dependent Ginzburg-Landau 
theory for multicomponent superconductors 


a transport time. This estimation does not depend on 
whether superconductor is in the clean or diffusive limit 
since a n is a normal state characteristic. The dimen¬ 
sionless conductivity in our units is given by <j n /<To ^ 
(vp / c) 2 (£f/T c )(Ht£f), where ep is the Fermi energy. 
Here the first factor is small since typically ( vp/c ) 2 ^ 
10 -5 in metals. The last two factors are large since 
(, Sf/T c ) ~ 10 2 and hence ( Htef ) ~ 10 2 (hrT c ). In this 
case in order to fulfil the clean limit conditions (tT c ) > 1 
we need to assume the values of dimensionless conduc¬ 
tivity CF n /(jQ > 0.1. 

In dimensionless units the TDGL equations read 


The time-dependent GL theory formally applies only 
in the gapless regime when either the inelastic electron- 
phonon relaxation r ~^ or spin flip rates are much 
larger than the superconducting pairing amplitude. We 
will consider here the former scenario which realizes in 
the vicinity of the critical temperature T c — T <C r^ 1 . 
In this case, the time-dependence is added to the generic 
Eqs.(A.15,A.17) in quite a standard way, following the 
derivation for a single-component case [50] 


[(Go + r — A _1 )A]^. = 

hix ( d 2 ie \ 


(B.l) 


+ |A,f A,-, 


Implementing the same reduction as in the stationary 
case, we obtain the system of two coupled TDGL equa¬ 
tions: 


^ ( d lie \ SB 47 r „ „ . 

Tk \Ft + ~n ip ) ^ k ~~W k ’ p VxB - a ™ E -^ 

(B.2) 

where k = 1, • • • , AT, T k = | Afc| 2 7 r/LBo/( 8 T c ) are the 
damping constants, cp is the scalar potential of electric 
field, a n is the normal state conductivity. The electric 
field is E = —c~ 1 d t A — V(/? where <p is electrostatic po¬ 
tential and for the superfluid current we have an expres¬ 
sion j s = —cSB/SA. In the following since |Ai 5 2| 2 ~ 1 
we put ri ,2 = T = 7 r hBft /(8T C ). 

In order to perform numerical simulations we normal¬ 
ize lengths by £0 = b,vp/T c , where vp is the average 
value of Fermi velocity, magnetic field by Bq = T c \Jvo/ p, 
free energy density by Fq = Bq, current density by 
jo = cBo/io, time by to = F/Hq, electric field by 
E 0 = £ 0 B 0 /(t 0 c), conductivity by a 0 = c 2 r/(^B$). In 
such units the electron charge is replaced by an effec¬ 
tive coupling constant e = irB 0 ^o/^o which parametrizes 
the regimes of extremely type-II and type-I supercon¬ 
ductivity at e C 1 and e 1 respectively. To es¬ 
timate the characteristic relaxation time we note that 
to = irh/8T c ~ 1 ps provided T c ~ 1 meV which is about 
10 K. 

The normal state electric conductivity a n in (B.2) 
can be defined from the Drude model. Up to unim¬ 
portant numerical coefficient resulting from the Fermi 
surface anisotropy, it yields a n ~ e 2 voVpT, where r is 


^Jf + 2 ej^ip k = -F- , VxB - a n E = j s (B.3) 

where E = —d t A — Vp and j s = —SB/SA. The free 
energy is given by 

E l n ^i1 2 + a Mo\ 2 + y IV’jl 4 } 

3 = 1 

+ fci2,a((n a ^i)*n a ^ 2 + C.C.) + 7 |^i| 2 IV’2| 2 

+ + c - c ) • (B.4) 

where II = V — 2 ieA. 


Appendix C: Numerical methods and starting guess 

We consider here the problem (B. 2 ) defined on a 
two-dimensional bounded domain D C R 2 with dfl its 
boundary. This problem is supplemented by the insu¬ 
lator boundary conditions, that no supercurrent flows 
through the boundary. The absence of supercurrent flow¬ 
ing through the boundary reads as n • II jj a = 0 with n 
the normal vector to dfl. The problem are then dis¬ 
cretized using a finite element formulation provided by 
the Freefem++ library [54]. Discretization of the inte¬ 
gration domain U is done using a (homogeneous) trian¬ 
gulation over D, based on Delaunay-Voronoi algorithm. 
The fields are decomposed on a continuous piecewise 
quadratic basis on each triangle. 


Numerical procedure 

In our simulations, we chose the relaxation factor for 
the components of the order parameter T = 1 and the 
normal state conductivity cr n = 0.1. We chose the gauge 
coupling constant e = 0.2 and the coefficients of the 
coupling matrix (A.9), that determine the coefficients of 
the Ginzburg-Landau functional (B.4), are p = 5 and 
A = 4.5. The coefficients in GL functional (B.4) are then 
consistently obtained using the previously defined rela¬ 
tions and the obtained values are f3i = 2, = 1.1079, 




S = 0.4645 and 7 = 0.929. Given the coefficients of gradi¬ 
ent terms in different bands are K W = 0.5, K^ = 0.05, 
K = 0.25 and using the relation (A. 17), the coefficients 
of the kinetic terms read as k\ = 0.55, £7 = 0.375 and 
ki 2 ,x = 0.217, while £ 72 , y = £7 2, x for s + is states and 
£72, y = —ki 2 , x for s + id. We are interested here in the 
effect of inhomogeneities of the order parameter induced 
by a local heating. In first approximation, temperature 
dependence within the Ginzburg-Landau theory is mod¬ 
elled by a linear dependence of the quadratic couplings 
on temperature: (T) = a^~ r, where r = (1 — T/T c ) is 

the reduced temperature and a&, a positive characteristic 
constants determined from the microscopic calculations, 
see Appendix A) . As the Ginzburg-Landau expansion can 
be justified only in the vicinity of T c , only a small interval 
of r can be investigated. 

The spatial modulation of the temperature is modelled 
as follow. Let To denote the outer boundary of the do¬ 
main (the sample’s boundary) and Ti an inner boundary 
(the heat source). In our simulations we considered var¬ 
ious shapes of Ti, but the results reported here are for 
T 1 to be a small circle. The temperature profile T(x) is 
determined by solving the stationary heat equation with 
a small damping factor fi 

A T(x) = fiT(x) and T(T^)=cte. (C.l) 

In our simulations, we set T(T 0 ) = To = 0.7T C for the 
sample’s boundary. On the other hand, the temperature 
of the source T s varies with time as: 

T -(')- r °+77 (1-771)) <a2) 

where denotes the floor function, T\ = 0.95T C is the 
maximal temperature of the source and tit is the num¬ 
ber of heat source temperatures. In our simulations, we 
chose t = 30. Fig. 5 shows a typical solution for the sta¬ 
tionary heat equation (C.l). Note that our choice here 
for considering the temperature profile given by the sta¬ 
tionary heat equation implies that we assume than heat 
transport occurs on time scales much smaller than the 
other time scales of the problem. 



0.95 

0.9 

0.85 

T/T c 

0.8 


0.75 


0.7 


The heat equation (C.l) is thus solved for a given tem¬ 
perature T s of the heat source (T 1 ). This define the tem¬ 
perature profile T(x,t) and thus the spatial modulation 
quadratic couplings c^(T) = di — r. The time dependent 
Ginzburg-Landau equations (B.2) are discretized using 
Euler’s forward method and iterated for a given inter¬ 
val At = NSt. The time step St depends on the various 
time scales of the problem. Here we chose St = 10 -2 and 
N = 8 x 10 3 . After A t , the procedure is repeated for a 
new temperature T s of the heat source given by (C.2). 


Starting guess for domain-walls 

The starting solutions are stationary solutions of the 
GL functional, for constant temperature. That is with 
no heat source. They can simply be ground state but 
also topologically non-trivial solutions in equilibrium. 
Domain-walls are topological excitations that are associ¬ 
ated with the spontaneous breakdown of a discrete sym¬ 
metry. These are field configurations that interpolate be¬ 
tween inequivalent ground-states that are disconnected. 
Here, we are in particular interested in domain-walls 
that interpolate between regions with inequivalent rela¬ 
tive phases between the condensates. When the ground- 
state breaks time-reversal symmetry, its complex conju¬ 
gate is not a gauge equivalent. That is, there exist no 
real number xo such that T^ = e ZXo To. If no such trans¬ 
formation exist, then T^ ^ T 0 and the configurations 
with ground state phases (p a and —ip a are disconnected 
and degenerate in energy. Domain-walls that interpolate 
between To and T^j are topologically protected as their 
unwinding would require to overcome an infinite energy 
barrier, see for example textbook discussion in [65, 66 ]. 
The domain-wall that interpolates between T^ and To 
can thus be parametrized by: 


V #*” 0 * U a exp 


iip a tanh 


x± - x 0 

Jdw) 

k Sa , 


(C.3) 


where Xq is the curvilinear abscissa that gives the posi¬ 
tion of the domain-wall, and x± is the coordinate per¬ 
pendicular to the domain-wall and £( dw ) determines the 
width of the domain-wall. 


Figure 5: (Color online) - Temperature profile obtained by 
solving the stationary heat equation (C.l) for a local heating 
of the sample. The inhomogeneous temperature profile is in¬ 
duced by the ring-like heat source shown show by the (purple) 
circle. 
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Appendix D: Description Additional Movie Material 

• Movie Anim-magnetic.avi: 

shows the evolution of the magnetic response that 
originates in time-varying local heating of the sam¬ 
ple. Detailed description is given in the caption of 
Fig. 6. 

• Movie Anim-imbalance.avi: 

shows the evolution of the electric response that 
originates in non-stationary local of the sample. 
Detailed description is given in the caption of 
Fig. 7. 
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Figure 6: (Color online) - Evolution of the magnetic response that originates in time-varying local heating of the sample. 
The overall time sequence runs for t/to = 2400. The leftmost panel shows the time evolution of the source’s temperature given 
by (C.2) and the red dot denotes the position in the time series. This corresponds to Figure 2 of the main text. The color 
plot shows the magnitude of the out-of-plane induced magnetic field B z in unit of the second critical field. The two upper 
panels show the response of an s + is superconducting state respectively with a domain wall and the homogeneous state. The 
two lower panels show the response of an s + id superconducting state respectively with a domain wall and the homogeneous 
(anisotropic) state. 
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Figure 7: (Color online) - Evolution of the electric response that originates in time-varying local heating of the sample. The 
overall time sequence runs for t/to = 2400. The leftmost panel shows the time evolution of the source’s temperature given by 
(C.2) and the red dot denotes the position in the time series. This corresponds to Figure 3 of the main text. The color plot 
shows the voltage induced in the normal detector by the charge imbalance due to the non-stationary heating. The two upper 
panels show the response of an s + is superconducting state respectively with a domain wall and the homogeneous state. The 
two lower panels show the response of an s + id superconducting state respectively with a domain wall and the homogeneous 
(anisotropic) state. 
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